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ABSTRACT We investigate the permeability of lipid membranes for fluorescence dyes and ions. We find that perme- 
ability reaches a maximum close to the chain melting transition of the membranes. Close to transitions, fluctuations in 
area and compressibility are high, leading to an increased likelihood of spontaneous lipid pore formation. Fluorescence 
Correlation Spectroscopy (FCS) reveals the permeability for rhodamine dyes across 100 nm vesicles. Using FCS, we find 
that the permeability of vesicle membranes for fluorescence dyes is within error proportional to the excess heat capacity. 
To estimate defect size we measure the conductance of solvent-free planar lipid bilayer. Microscopically, we show that 
permeation events appear as quantized current events. Further, we demonstrate that anesthetics lead to a change in 
membrane permeability that can be predicted from their effect on heat capacity profiles. Depending on temperature, the 
permeability can be enhanced or reduced. We demonstrate that anesthetics decrease channel conductance and ultimately 
lead to 'blocking' of the lipid pores in experiments performed at or above the chain melting transition. Our data suggest 
that the macroscopic increase in permeability close to transitions and microscopic lipid channel formation are the same 
physical process. 
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Introduction 

Lipid membranes are regarded as insulators that are practically 
impermeable to ions and larger charged molecules. This is cru- 
cial for the Hodgkin-Huxley model for the action potentiall 1 1 
but also for the interpretation of many microscopic permeation 
events associated with proteins. In the vicinity of the melt- 
ing transition, however, membranes become permeable to wa- 
ter |2), ions OHH and even to large molecules. This may be 
of some relevance since many biological membranes are in fact 
found close to such a regime J5]|6l. It has been known since 
the early 1970s that the permeability of lipid membranes ap- 
proaches a maximum in the chain melting regime of lipids in 
the absence of proteins |7], and this is also likely to be so for 
biological membranes in the transition regime. The reason for 
this enhanced permeability is the enhanced area fluctuations 
that lead to a maximum in the lateral compressibility (8][9]. 
Nagle and Scott proposed that the changes in lateral com- 
pressibility lead to a facilitation of pore formation since the in- 
creased compressibility lowers the work necessary to create a 
membrane defect. In the melting transition one also finds do- 



main formation. Papahadjopoulos and coworkers Q suggested 
in their seminal paper that permeation events take place at do- 
main boundaries. This view has also been adopted by Mourit- 
sen and coworkers ifTOlfTTI . At domain boundaries, fluctuations 
in lipid state are maximal, and this suggests that domain inter- 
faces are especially leaky. However, this view is not identical 
to that of Nagle and Scott [8 ] since the overall length of domain 
interfaces is not identical to large fluctuations. Unfortunately, 
the overall quality of permeation data for membranes so far has 
been low. This is partially due to difficulties in maintaining 
constant temperature during optical experiments. Melting tran- 
sitions of single lipid LUV display half widths of less than one 
degree. For a quality permeability experiment on such a sys- 
tem, an accuracy of at least 0.1 K is required. 
While the measurements mentioned above refer to ensembles, 
i.e. dispersions of vesicles, one can also monitor permeability 
for ions on microscopic scales. This is done either by using 
black lipid membrane (BLM) spanning a small hole in a teflon 
film with a diameter of 50-100 /im, (Montal-Muller technique, 
JT2]) or in patch clamp measurements where one uses mem- 
brane patches of the order of 1 fim 2 . Surprisingly, measure- 
ments of the conductance of pure synthetic lipid membranes 
can display quantized steps very similar to those reported for 
ion channel proteins. While there is, to our knowledge, no con- 
vincing explanation of why such currents should be quantized, 
there are several publications documenting such channel-like 
behaviour. The first paper on this is from Yafuso et al. Ifl3l 
from 1974 on oxidized cholesterol films, which appeared two 
years before the famous paper by Neher and Sakmann lfT4l on 
quantized currents through the acetylcholine receptor. Antonov 
and collaborators 1 3j were the first to show that these lipid mem- 
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brane conduction events are strongly enhanced in the phase 
transition regime. Using pure synthetic lipids, they could switch 
permeation on and off merely by shifting the temperature by a 
few degrees. This basically rules out the possibility that one 
is observing proteins or peptide impurities. Kaufmann lfl5l 
[TBI [171 showed that the occurrence of such lipid ion channels 
depends not only on temperature but also on other intensive 
variables such as lateral tension and pH (i.e., the chemical po- 
tential of protons). Along the same lines, Antonov and col- 
laborators 1T81 as well as Gogelein & Koepsell lfl9l showed 
that channel-like events are influenced by calcium concentra- 
tion (i.e., the chemical potential of calcium). Furthermore, they 
depend on voltage fl20| . All of these changes are also known 
to influence transition temperatures in lipid membranes. Sev- 
eral other publications describe quantized ion currents in pure 
lipid membranes in the absence of proteins |2"T1 122"1 |23l Pfl . In 
spite of the putative relevance of such findings that are well 
documented and which show the coupling to thermodynamical 
variables convincingly, this phenomenon is not well known to 
a broad community. 

It is well-established that anesthetics change the melting be- 
havior of membranes. In particular, they lead to a lowering 
of the transition temperature 11241 l25l 1261 . Thus, it seems ob- 
vious to study the effect of anesthetics on permeability. It is 
an interesting historical note that Charles Ernest Overton, who 
first demonstrated the general action of anesthetics and its re- 
lationship to membrane solubility (Meyer-Overton rule, |27|) 
was also the first to studied the permeability of membranes to 
anesthetics (Overtonl895). 

In this publication we compare ensemble permeation experi- 
ments on lipid vesicles with ion conductance measurements in 
planar lipid membranes. We use the FCS technique to measure 
the ensemble permeation rate for fluorescent dyes. In particu- 
lar, we study the influence of temperature and the presence of 
anesthetics on permeability. Then, we determine the permeabil- 
ity for ions in a BLM setup. Here, we vary temperature, voltage 
and anesthetic concentration. The goal is to demonstrate the 
thermodynamic couplings that lead to changes in permeability. 
In particular, we show that anesthetics influence the permeabil- 
ity in a manner closely related to their influence on the melting 
of lipid membranes that can be explained by freezing-point de- 
pression. 

Materials and Methods 

Chemicals Decane and 1-octanol were purchased from Fluka 
Chemie AG (Deisenhofen, Germany), N-hexadecane, chloro- 
form and methanol were obtained from Merck (Hohenbrunn, 
Germany), n-pentane was provided by BDH (Poole, UK) and 
potassium chloride from J.T Baker Analyzed (Deventer, Hol- 
land). l,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 
l,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG) and 1,2- 
dioleoyl-sn-glycero-3- phosphocholine (DOPC) were purcha- 



sed from Avanti Polar Lipids (Birmingham, AL) and used with- 
out further purification. Rhodamine 6G and tetramethylrho- 
damine dextrane were purchased from Invitrogen/Molecular - 
Probes (Carlsbad, Ca). For all experiments, MilliQ water (18.1 
Mfl) was used. 

Fluorescence correlation spectroscopy: 
Sample preparation. We prepared unilamellar lipid vesicles 
(LUV) from DPPC:DPPG=95:5 mol/mol with an Avestin ex- 
truder (Avestin Europe, Mannheim, Germany) using a filter 
with a pore size of 100 nm diameter. Extrusion was performed 
in the presence of high concentrations of Rhodamine 6G flu- 
orescence markers (50 /iM). Subsequently, remaining free dye 
was removed on a G50 Sephadex column at temperatures below 
the melting temperature of the lipids where the membranes are 
nearly impermeable. We used vesicles made of a DPPCDPPG- 
=95:5 mixture. DPPG is negatively charged and was added to 
prevent the aggregation of vesicles. All experiments were per- 
formed in a 200mM NaCl buffer. At this salt concentration, 
vesicles do not aggregate but R6G (positively charged) does 
nevertheless not associate with the membranes. 
Fluorescence correlation spectrometer: Our inverted micro- 
scope setup has been describe in detail in l28l . We used a lin- 
early polarized continuous wave 532 nm Nd:Yag laser (Laser 
2000, Wessling, Germany) with a power of 5mW. Further, we 
used a 1 .20 NA 60 x water immersion objective (Olympus; UP- 
LAPO) and a confocal setup with pinhole sizes of 30 /im. The 
fluorescence signal was detected by a SPCM-AQR-13 avalanche 
photo diode (Laser Components, Olching, Germany). Time 
scales were calibrated with a Rhodamine 6G solution at 296 
K with a known diffusion coefficient of D = 3 • 10~ 6 cm 2 /s at 
22°C. The signal from the APD was analyzed using a FLEX- 
5000/fast correlator card by Correlator.com (Bridgewater, NJ). 
In order to adjust temperatures in the FCS experiment, both 
sample and microscope objective were temperature-controlled 
using a HAAKE DC30 K20 (Waltham, USA) waterbath. Ad- 
ditionally, the entire setup was heated by a ventilating radiator 
so that ambient temperature was close to the experimental tem- 
perature at the objective. 

Calorimetry: DSC experiments were performed on a VP- 
DSC Calorimeter (MicroCal, Nort-hampton, MA, USA) with a 
scan rate of 5°C/h. The calorimetric experiment related to the 
FCS experiment (Fig. [2]) was performed on the same sample 
as the one used in the FCS experiments. For the DCS exper- 
iments related to the BLM experiments (Figs. |4] and |5), the 
lipid samples were prepared by pre-dissolving in chloroform 
and drying the solvent under vacuum overnight. The dried lipid 
mixtures were dispersed in MilliQ water to a final concentra- 
tion of 20mM. The buffer used was the same as for the BLM 
experiments. Before filling the calorimeter, the solutions were 
degassed for 10 minutes. 

Black lipid membranes: Planar bilayers were formed over a 
round aperture in a Teflon film of 25 /i m thickness, dividing two 
compartments of a teflon chamber embedded in a brass block 
that could be heated by a circulating water bath. The aperture 
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of «80 /im radius was punctured by a steel needle (experiments 
in Figs. [5] and |6]l or by an electric spark (Fig. @). Lipid mix- 
tures were prepared from chloroform-methanol solutions. The 
samples were dried under a weak flow of nitrogen gas/air and 
placed under vacuum overnight to remove the residual solvent. 
The aperture in the Teflon film was prepainted with 5% hexade- 
cane in pentane. The BLMs were painted with DOPGDPPC 
2:1 lipid solutions in decane/chloroform/methanol 7:2:1 and 
formed following the method described by Montal and Miiller 
lfl2l . The two compartments of the Teflon block were filled 
with unbuffered 150mM KC1 (pH 6.5). Lipid solution (25mg/ml) 
was spread on the buffer surface in each compartment (approx. 
3 /xl on each side). Ag/AgCl electrodes were placed into both 
compartments of the chamber. After 15-30 minutes to allow 
for the evaporation of the solvent, the water level of the com- 
partments was lowered and raised several times until a bilayer 
was formed over the hole. The formation of a BLMs was con- 
trolled visually and by capacitance measurements (with trian- 
gular lOOmV voltage input pulse). The specific capacitance of 
the membranes was found to be 0.9- 1 .2 /iF/cm 2 . For the exper- 
iments with octanol, 15% v/v octanol in methanol solution was 
prepared and added symmetrically on both sides of the Teflon 
chamber. 

Conductance measurements were performed on an Axopatch 
200B amplifier in voltage clamp mode connected to a DigiData 
1200 digitizer (both Molecular Devices, Sunnyvale, CA, USA). 
Current traces were filtered with 1kHz low -pass Bessel filter 
and recorded with Clampex 9.2 software (Axon Instruments) 
on the hard drive of the computer using an AD converter with 
a time resolution of 0. 1ms. The data was further analyzed with 
Clampfit 9.2 and low-pass filtered with Bessel (8-pole) filter at 
a cut-off frequency of 300Hz. Temperature was controlled by a 
HAAKE DC30 K20 (Waltham, USA) waterbath and a thermo- 
couple (WSE, Thermocoax). 

The BLM experiments in Fig. |4]were performed in the labora- 
tory of M. Winterhalter in Bremen, the experiments shown in 
Figs. |5]and|6]were performed using a very similar setup in the 
laboratory of T. Heimburg in Copenhagen. 

The FCS correlation function of vesicles with variable con- 
tent of fluorescence dyes: The correlation function of a freely 
diffusing dye in a microscope focus of Gaussian cross section 
is given by: 



correlation as a square 
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where the brightness of particle species i is given by B L = 
KitiQi l29l . Here, Ki is the efficiency of the fluorescence de- 
tector, 6j is the molar extinction coefficient of the fluorophore 
at the wavelength of excitation, and Qi is the quantum yield of 
the fluorophore. For fluorophores of equal brightness, the Bi& 
cancel. In the presence of a single dye species, one obtains eq. 

m 

Permeation through vesicles studied by FCS In this work we 
wish to monitor the permeation of fluorescence markers from 
the inside of a lOOnm vesicle to the outside. These vesicles 
are smaller than the microscope focus. Even though all mark- 
ers have the same brightness in this particular experiment, a 
vesicle that contains n dye molecules is n time brighter than 
a single dye molecule (neglecting quenching effects). In par- 
ticular, vesicles do not all contain precisely the same number 
of fluorophores. The dye entrapped in vesicles and that free 
in solution is rhodamine 6G. Thus, the difference in brightness 
between freely diffusing single dye molecules and the vesicles 
is directly proportional to the number of dye molecules trapped 
within the vesicle. Therefore, in the following only the number 
of dye molecules per vesicles will appear. Not all vesicles will 
contain precisely the same number of molecules. From statis- 
tics, one expects a Poisson distribution of the dye content of the 
vesicles: _ _ 
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where B is the mean number of dyes per vesicle. If one inserts 
this distribution into eq. [2] one obtains 
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If one has several markers of similar 'brightness', the correla- 
tion functions add, and the number of diffusing particles can be 
deduced from the relative amplitude of the steps in the correla- 
tion function. The situation becomes more complicated if dif- 
ferent objects with different brightnesses are present or if there 
is a distribution of brightnesses since the brightness enters the 
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The detailed derivation can be found in the M.Sc. thesis of A. 
Blicher 1301 . This equation contains five unknown parame- 
ters that must be determined from experiment. The correlation 
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Fi gure 1 : Principle of the measurement of permeation of the fluorescence dye R6G in the FCS experiment. Left: Schematic 
drawing. Dye molecules trapped in vesicles smaller then microscopic focus size diffuse with the characteristic time constant of 
the vesicles. With time, molecules permeate through the vesicle membranes and diffuse with the characteristic time constant of 
free molecules. Center: The correlation profiles of this situation show two correlation steps corresponding to free molecules 
and vesicles. The relative step size changes as a function of time such that the number of free molecules increases and the 
number trapped in vesicles decreases. Five experiments are shown that were performed after the times indicated. Right: The 
mean number of dye molecules per vesicle (calculated from the mean brightness, see text) decreases exponentially with time. 
Data shown have been recorded on DPPC.DPPG (95:5 mol/mol) LUV in 200mM NaCI at temperatures of 39.5° C (center 
panel) and 36.6° C (right panel) below the chain melting transition maximum (see Fig. |2j- At 36.6°C the time constant of 
permeation is 4200s. Each data point corresponds to an FCS experiment as shown in the center panel. 



function yields the two diffusion correlation times trqg an d 
Ty of Rhodamine 6G and the vesicles, respectively. The mean 
number of free dye molecules in the focus, (Nrqg), the mean 
number of vesicles in the focus, (Ny), and the mean number of 
dye molecules per vesicle, B, can be obtained from the absolute 
fluorescence intensity of the uncorrected data, and the magni- 
tude of the two steps in the correlation function (see FigJT}. In a 
time-dependent experiment, the mean number of dye molecules 
per vesicle will decrease. The change of B with time is propor- 
tional to the permeation rate through the vesicular membrane. 

Theory 

The connection between compressibility, heat ca- 
pacity and permeability 

The permeability of the membrane is maximal in the transition 
regime. This phenomenon is related to the well-known phe- 
nomenon of the opalescence of fluids close to critical transi- 
tions. In this regime, fluctuations in density may approach the 
length scale of visible light. Einstein [31 J calculated the work 
required to generate density fluctuations on the length scale of 
light. In particular, he determined the work necessary to change 
the density of box of size L 3 by moving volume into another 
box of equal dimensions. Close to transitions this is consider- 
ably easier because the volume compressibility becomes high. 



In fact, one can see the increased fluctuations in liquids close 
to their transition in optical experiments. The problem of per- 
meability changes and of pore formation is of a similar nature 
since it requires moving lipids into another region of the mem- 
brane in order to create a pore. In particular, the likelihood 
to form pores or holes in the membrane is closely related to 
fluctuations in density. The problem can now be posed as fol- 
lows: How much work must be performed to move a lipid from 
a membrane segment of dimensions L 2 into another region of 
similar size? 

The work necessary to compress a membrane (i.e., the Helmholtz 
free energy change) by an area a is a quadratic function of the 
area that has the form 

AWr(fl) = ^(i) ' A ° (5) 

where Aq is the total area of the membrane, a is the size of 
a defect and is the isothermal area compression modulus 
(e.g., [8]). The compressibility of a membrane is given by 

A _ 1 (dA\ 1 

Kt ~ To\m) T -K* (6) 

where is the isothermal area compressibility. The compress- 
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ibility is also proportional to the area fluctuations |]9| 



Permeability of LUV measured by FCS 
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while the heat capacity is proportional to the fluctuations in en- 
thalpy (9): 
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The empirical finding that, close to transitions, the changes in 
area and the changes in enthalpy are proportional functions of 
temperature, i.e., AA(T) = ~/ A ■ AH{T) |9]|32][33]|33, leads 
to 
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Nagle and Scott [8 | assumed that the permeability can be ex- 
pressed as a series expansion with respect to the area fluctu- 
ations that are also proportional to the lateral compressibility. 
The resulting permeability obtained by Nagle and Scott is 
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higher order terms 
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where ciq and (12 are constants to be determined from experi- 
ment. Using expressionfTOlvields 
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where cq and C2 are also constants to be determined from ex- 
periment. This is a simple expression directly coupled to the 
calorimetric experiment. In particular, it predicts that the change 
in permeability is proportional to the excess heat capacity, which 
is easily measured. This approach relies on area fluctuations 
and does not involve a molecular model for a pore. In partic- 
ular, we do not consider a line tension around the pore. The 
experimental results below will justify this approach. 
The results of this section also immediately imply that every- 
thing that changes heat capacity profiles will also change the 
permeability in a coherent and predictable manner. The similar- 
ity of this problem to critical opalescence implies that it should 
be possible to detect the critical nature of pore formation in op- 
tical experiments. 

Results 

We have investigated two different aspects of lipid membrane 
permeability: the macroscopic permeability of lipid vesicles for 
fluorescence markers and the permeability of planar lipid mem- 
branes for ions. While in the first setup one obtains an ensemble 
average over many permeation events, in the second experiment 
the permeation events allow us to draw conclusions regarding 
the possible size of the defects. In particular, we investigate the 
correlation with anesthetics. We show that such drugs influence 
the permeability in a coherent and predictable manner. 



Fluorescence correlation spectroscopy is a single molecule tech- 
nique in which the diffusion constant and particle concentra- 
tion can easily be determined on the length scale of the micro- 
scope focus. Rhodamine 6G fluorescence dyes were entrapped 
in lipid vesicles of 100 nm diameter. Both, vesicles and dyes 
are smaller than the microscope focus diameter. However, the 
diffusion constant of the singular dyes is much higher than that 
of the much larger vesicles. While the typical dwell time of a 
single dye molecule in our FCS setup is on the order of 100 
/is, the dwell time of a fluorescing vesicle is of order 10ms. If 
both single dye molecules and fluorescent vesicles coexist in 
a solution, the correlation function will show two steps with 
amplitudes related to the concentrations of dyes in vesicles and 
free dyes. 

Initially, all R6G molecules are trapped in the vesicles. With 
time, they leak out of the vesicles as shown schematically in 
Fig- HI (left). This process can be followed in repeated FCS 
experiments (Fig. Q] center), and the relative fraction of dye in 
vesicles and as free molecules can be deduced from the ampli- 
tudes of the steps in the correlation function following the pro- 
cedure given in eq. |4]in the Materials and Methods section. As 
an example, we show the leakage of R6G out of vesicles made 
of 95 mol% DPPC and 5 mol% DPPG measured at 39.5°C, 
which is about 2 K below the heat capacity maximum of the 
membrane melting profile. Under these conditions permeation 
of dyes through the vesicular walls is quite slow. One can see 
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Figure 2: Permeation rate of R6G through DPPC:DPPG 
(95:5) LUV (200mM NaCI) as a function of temperature com- 
pared to the heat capacity profile of the identical sample. The 
permeation rate is closely related to the heat capacity as pre- 
dicted by eq. \12\ 
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Figure 3: Effect of the anesthetic octanol on the permeation rate. Left: Permeation of R6G through DPPC.DPPG (95:5) 
LUV at 37.1° C after the sudden addition of 1 -octanol at the time indicated by the arrow. Permeation increases dramatically. 
Center: Permeation in the absence and the presence of 3.5mol% and 7 mol% octanol in the membrane at 37.6° C. Right: Heat 
capacity profiles of the lipid LUV in the absence of octanol (cf. Fig. |2j and in the presence of7mol% octanol. In the presence of 
octanol the melting events occur at lower temperature. The filled circles indicate the permeation rates at 37.6° C in the absence 
of octanol and in the presence of 7 mol% octanol. The increase in permeability caused by octanol correlates nicely with the 
increase in heat capacity at 37.6° C. 



two correlation steps. While the correlation times remain con- 
stant, the relative fraction of slow and fast diffusion species is 
changing. FigQ](right) shows that the progress of leakage for a 
sample at 36.6°C (« 5 K below the transition maximum) yields 
a single exponential decay with a time constant of 4800 seconds 
and a permeation rate constant of 2.1 • 10~ 4 s _1 . 

We recorded the permeation rate constants as a function of tem- 
perature and compared them to the heat capacity profile (Fig. 
[2]). It was found that the permeation rate of R6G was, within 
experimental accuracy, linearly related to the excess heat ca- 
pacity as predicted by eq. [12] In the temperature range from 
36.5° to 41.4°C, the permeability changed by a factor of 50. 
Compared to temperatures below the pretransition temperature 
(< 34°C), the change is several orders of magnitude (data not 
shown). Close to the heat capacity maximum the permeation 
rate is about 0.01 s _1 , which corresponds to a time constant 
of 100 seconds. The dwell time of a vesicle in the microscope 
focus is of order 10 ms. To obtain a correlation function of 
reasonable quality, one needs to record a fluorescence noise 
trace of at least 10 s. Thus, the permeation rate at heat capacity 
maximum is close to the resolution limit of our method. Only 
two data points have been measured above the melting regime. 
These points are difficult to obtain because, during equilibra- 
tion of the sample to the experimental temperature, it passes 
through the phase transition regime and most of the dye leaks 
out prior to the start of the experiment. The permeation rates 
for temperatures above 41.4 C given in Fig. [2]represent the 
single exponential decay of the remaining fluorescence mark- 
ers (of order 10% of the original dye concentration). This slow 
decay is absent in the data at the transition point. We conclude 
that the permeability returns to values similar to those below the 
transition range, in agreement with previous studies. A control 



experiment with tetramethylrhodamine dextrane (MWw 3000 
g) at the melting temperature of 41.4 °C showed that this large 
dye remained inside of the vesicles on time scales of several 
hours (data not shown). This indicates that the vesicles stay in- 
tact when passing through the melting regime and do not rup- 
ture. 

Influence of anesthetics: In 2007 we showed that anesthet- 
ics act as freezing-point depressants l26l . At critical anesthetic 
dose (2.6 mol% of anesthetics in the membrane), one obtains 
a melting point depression of about 0.6 degrees. Since anes- 
thetics shift the heat capacity profile by a known value, one 
expects a predictable change of anesthetics on permeability. In 
particular, if one measures permeability below the melting tem- 
perature of the vesicles, anesthetics should increase the perme- 
ability since the melting events are moved towards experimen- 
tal temperatures. If one repeats the same experiment above the 
melting temperature, addition of anesthetics should decrease 
the permeability since the transition events are moved away 
from experimental temperatures. This was investigated using 
the anesthetic 1 -octanol. Fig [3] (right) shows the effect of the 
addition of 3.5 mol% percent of octanol (calculated from the 
total octanol content and a lipid/water partition coefficient of 
200 1 35 1) on the heat capacity profile. It shifts by about 0.9 
K to lower temperatures, in agreement with the freezing-point 
depression relation derived in 11261 . Upon addition of octanol 
to the vesicle dispersion at 37.6 °C, the permeation rate in- 
creases dramatically (Fig(3j left. The time of octanol addition 
is indicated by an arrow). Dye release for different amounts 
of octanol is shown in (Figj3] center). The permeation rates 
are 0.24 • 10~ 3 s~ 1 in the absence of octanol, 1.19 • 10 _3 s _1 
in the presence of 3.5 mol% octanol and 19.8 • 10~ 3 s _1 in the 
presence of 7 mol% octanol. In Fig[3] (right), the permeation 
rate with and without octanol are compared with the heat ca- 
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Figure 4: Spontaneous quantized currents through BLMs made of a DOPC:DPPC=3:l mixture (150mM KCl) as a function 
of temperature. Left: Heat capacity traces with a heat capacity maximum around 17° C. At 30° C is above the melting regime 
in its fluid phase. Center: Current traces at 17° C (40mV) and 30° C (150mV). One only finds current steps in the lipid melting 
regime. Right: Current histograms for the two temperatures. At 17° C two maxima in the histogram can be seen. 



pacity changes. At 37.6°C the permeability increase caused by 
octanol is about 5-fold and matches exactly the change of the 
excess heat capacity caused by octanol at the same temperature. 



Permeation of ions through black lipid membranes 

It has been shown previously J3] H) that permeation of ions 
through BLMs is most prominent in the phase transitions regime 
Further, conductance changes occur in quantized steps simi- 
lar to that reported for channel forming proteins or peptides 
H3] [3] HU SI [H] HSl HU [I71. In order to stay closer to room 
temperature, we used lipid mixtures in the BLM experiments 
which differed from those of the FCS experiments. BLMs be- 
low the transition are typically not stable. In Fig|4]we show data 
recorded on a DOPC:DPPC=3:l (mol/mol) mixture (150mM 
KCl, ImM Tris, pH 7.4). This mixture displays a broad melt- 
ing profile with a maximum at 17°C (Fig. [4] left panel). Above 
22°C the mixture is in the fluid phase. We compared current 
measurements at 17°C and 30°C (Fig. |4] center panel). At 
17°C and a transmembrane voltage of 40m V, one sees discrete 
current steps with an amplitude of about 35pA (see histogram 
in the right hand panel) and typical channel open times of sev- 
eral 10ms. At 30°C one does not see any current steps even 
at higher voltages of 150mV. This supports the notion that the 
transition is the regime in which such events can be found. Our 
observation of well-defined conductance steps indicates that 
this particular lipid composition is able to form well-defined 
defects even in the absence of peptides or proteins (e.g, |fT9l ). 
Fig. |5](top) shows traces from a DOPC:DPPC=2:l (150mM 
KCl, 60mV, pH « 6.5, 19°C). This mixture has a heat ca- 
pacity maximum close to 19°C. The figure shows that one can 
find traces with at least four equally spaced conductance steps, 
probably corresponding to several lipid channels opening si- 
multaneously. Such multiple events seem to be more frequent 



when one is closer to the transition regime or if the voltage is 
higher. The bottom panel of Fig. shows a current-voltage 
relationship for the membranes used in the top panel recorded 
at 19°C. It represents the average current for all events where 
I m is larger than zero. It is linear in the regime between - 
150mV and +150 mV. If |V| >150 mV, the overall conduc- 
tance increases, probably due to a significant contribution of 
multiple conductance steps induced by higher voltage. Note 
that the electrostatic potential is an intensive thermodynamic 
variable that influences the state of the membrane (see discus- 
sion). From the slope in the linear regime one can calculate 
a channel conductance of 300±80 pS. Assuming a pore in the 
lipid membrane filled with electrolyte J4|, one calculates a pore 
radius of about 0.75nm which is practically identical to the di- 
ameter of a single lipid. (DPPC has a cross-sectional area of 
0.63nm 2 in the fluid state J9) corresponding to a diameter of 
about 0.8nm.) This size is similar to the radii reported for pro- 
tein channels, e.g., 0.45-0.5 nm for the voltage-gated potassium 
channel Kvl.2 11361 or 0.3nm for the acetylcholine receptor 
[1371 . The fact that lipid ion channels display larger conduc- 
tance at high voltage could be considered as a voltage gating. 
However, it seems likely that this effect rather reflects the influ- 
ence of voltage on the transition. 

In the section above on dye permeation through vesicles it was 
shown that anesthetics (octanol) influence permeation due to 
its influence on the heat capacity and hence the lateral com- 
pressibility. In Fig. [6] (left panel) the heat capacity profiles of 
DOPC:DPPC=2: 1 membranes are shown in the absence and the 
presence of 1-octanol. For the DPPC:DPPG=95:1 membranes 
in Fig. [3] the melting events are shifted towards lower temper- 
atures. We performed BLM measurements at 19°C, which is 
close to the transition maximum. Thus, one expects that adding 
octanol moves the membranes out of their transition regime and 
according to eq. Q~2] the permeability should decrease. That 
this is in fact the case is shown in Fig. |6](center and right). At 



7 



< 



40- 
0- 



DOPC:DPPC 2:1, 150mM KCI 
pH~tf.5,60mV, 19°C 




voltage [raV] 

Figure 5: Currents and current-voltage relationship for BLMs 
made of a DOPC-DPPC=2:l mixture (150 mM KCI, 19° C, 
pH~ 6.5). The membrane is in its melting regime (cf. Fig. 
[6] left panel). Top: Current traces at a voltage of 60mV. One 
can see up to five current levels with equal distance between 
adjacent levels. Bottom: Current-voltage relation for the above 
membranes. The mean amplitude of current events is shown as 
a function of voltage. It is linear in the range between -150mV 
and +150mV. If \V\ is larger than 150 mV, the relationship is 
non-linear due in part to several simultaneous open events in 
the membrane. A total of 328 current traces of 30 second dura- 
tion was used for this analysis. 



210mV, the currents in the absence of octanol show two max- 
ima around 100 and 200pA. In the presence of 7.9 mol% oc- 
tanol the current steps range between 20 and 40pA, while in 
the presence of 15.9 mol% octanol the current amplitudes are 
on the order of 15pA. At higher octanol concentrations (23.8 
mol%), conductance events disappear completely. This reflects 
the decrease in heat capacity at 19°C induced by increasing 
octanol concentrations. Visual inspection of the current traces 
also seems to indicate that the mean open time per channel de- 
creases with increasing octanol content. 



Discussion 

In this paper we investigated the permeability of lipid mem- 
branes for dyes and ions using FCS and BLM measurements. 
In agreement with earlier studies HO HUB), we found that the 
permeability for both dyes and ions is maximum in the melting 
regime of the membranes. The FCS experiments showed that 
the permeation rate is, within experimental error, proportional 
to the heat capacity of a DPPC-DPPG=95:5 membrane. This 
can be understood by the large fluctuations in area near the tran- 
sition. This leads to a high compressibility and, consequently, 
the work required to create a pore by thermal fluctuation is 
small (see Theory section). Several degrees below the transi- 
tion the membrane was practically impermeable to dyes. Using 
the black lipid membrane (BLM) technique HI 21 . we measured 
ion currents through lipid pores. The findings are in agree- 
ment with the dye permeation experiments. We showed that 
one finds quantized currents in the melting regime which are 
absent at higher temperatures. Due to their similarity with pro- 
tein channel currents, one can consider these events as lipid ion 
channels. Leirer et al. 11381 showed that the mean conductance 
of a membrane of a lipid mixture indeed correlates well with 
the heat capacity (in their case of a DCisPCrDMPC mixture). 
Glaser et al. |39l suggested that at sufficiently high voltages 
stable hydrophilic pores of 0.6-lnm can exist. This is in agree- 
ment with our findings that suggest pore sizes on the order of 
one lipid cross-section. The ionic radius of K + is 0.138nm, the 
radius of CT~ is 0.181nm and the size of the dye molecule is on 
the order of lnm. In contrast to protein channels, the number 
of lipid channels is not fixed. They are merely fluctuating de- 
fects in membranes, and their occurrence depends on the state 
of the membrane as controlled by the intensive thermodynamic 
variables of the system (temperature, pressure, electrostatic po- 
tential, and the chemical potential of the components). We sug- 
gest that, if number and lifetime allow it, such lipid pores may 
fuse to form pores that are large enough to be permeable to 
large molecules. If one is precisely at the transition, this may 
already occur at zero or small voltages. It is known that at high 
voltages one finds the well-documented phenomenon of elec- 
troporation |40l . At and above 500m V, long-lived pores can be 
produced in biological cell membranes. This effect is used to 
facilitate the transport of drugs into cancer cells (4T|. Molec- 
ular dynamics simulations suggest that stable pores of nm-size 
can exist after applying transmembrane voltage pulses Il42ll43l . 
We have shown here that at voltages above 150mV the current- 
voltage relationship becomes non-linear so that currents are 
much larger. A similar finding on DCi5PC:DMPC mixtures 
was reported by (38). We further showed the existence of pores 
of a size sufficient to be permeable for small fluorescence dyes. 
One may well ask if the quantized currents found at lower volt- 
ages and electroporation found at high voltages are, in fact, the 
same phenomenon with the difference that the pore sizes are 
larger at high voltage. This matter requires further investiga- 
tions. In our experiments, the BLMs typically rupture above 
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Figure 6: Influence of the anesthetic octanol on membrane currents through BLMs made of a DOPC:DPPC=2:l mixture (150 
mM KCl, 19° C, p/f~ 6.5) at a voltage of 210 mV. Shown are four different octanol concentrations, mol% (bottom panels), 7.9 
mol% (2nd from bottom), 15.9 mol% (2nd from top), and 23.8 mol% (top panels) given as mol% of molecules in the membrane. 
Left: Heat capacity traces for three octanol concentrations (trace for 23.8mol% not shown). Increasing octanol leads to a shift 
of the profiles towards lower temperatures. This effect is known as freezing-point depression H26V . In the absence of octanol, 
the experimental temperature is at the c p -maximum, while it is above the c p -maximum in the presence of octanol. Center: 
Representative current traces at the four octanol concentrations. The current amplitudes and the total time in the open state 
become smaller with increasing octanol concentration. Right: Current histograms showing several peaks. At higher octanol 
concentration, the peak positions are at lower voltages. At 23.8mol% current events disappear completely. Note the different 
scales of the current axis in the center and the right panel. 



300m V. It should be noted here that one major disadvantage of 
the BLM technique is that the Teflon hole is prepainted with 5% 
hexadecane or decane in pentane. Thus, the membranes contain 
some residual alkanes that may influence the state of the mem- 
brane to a degree that may vary between different experiments. 
It would be of great interest to develop the Montal-Miiller tech- 
nique so that no additional organic solvent is required — or to 
repeat these studies with the patch clamp technique, as for ex- 
ample done in IfTTl . 

We have shown above that the presence of a melting transi- 
tion and the accompanying area fluctuations are necessary for 
permeation effects. Whatever changes the state potentially ef- 
fects the permeability. Here, we studied in particular the ef- 
fect of anesthetics on the permeation of dyes in FCS and ions 
in the BLM technique. We found that in both techniques the 
permeability changes in agreement with the changes in heat ca- 
pacity. The lipid ion channels showed smaller conductances in 
the presence of anesthetics. In a recent paper |26| we inves- 
tigated the effect of anesthetics on membranes that lower the 
melting point of membranes. Here, we have shown that, below 
the melting transition, the presence of anesthetics can increase 
the permeability because the transition temperature is shifted 
towards the experimental temperature (see Fig. [3]). If one is 
above the transition, the presence of anesthetics lowers the per- 
meability because the transition events are moved further away 
from experimental conditions (see Fig. [6j. In 1261 , we derived 
the following equation for the difference in free energy between 



fluid and gel phase of a membrane 
AG 



Arr , rm-r RT T 

±m tin J 77). 



+ f(l^) + f(pca^) + /(*) + ...) (13) 

where AH is the transition enthalpy, T m is the melting temper- 
ature, xa is the molar fraction of anesthetics in the membrane, 
Ap is the hydrostatic pressure and jv is a parameter describ- 
ing the relation between volume and enthalpy changes (jv = 
7.8 x 10 -10 m 2 /N). This equation contains one term describing 
the effect of temperature, a term for the influence of anesthet- 
ics and a further term describing the effect of hydrostatic pres- 
sure. For each change in an intensive thermodynamic variable, 
one can add further terms into this equation, for example for 
pH, calcium concentration or voltage, denoted here as J{fL H + ), 
f([iCa 2 +)> an d /O^)- We suggested in |26| that whenever the 
value of the potential AG is the same, the anesthetic action 
is also the same. Similarly, one can state here that whenever 
the value of AG is the same, one expects the same permeabil- 
ity. Changes in different variables may compensate, for exam- 
ple that of pressure and the anesthetics concentration known as 
pressure reversal of anesthesia ll26l . Therefore, it is to be ex- 
pected that the decrease in permeability observed in our BLM 
measurements would be compensated by hydrostatic pressure. 
In biological preparations that are found slightly above the tran- 
sition of their membranes, one would expect an increase in con- 
ductance of lipid ion channels with increased pressure. In this 
context, it is interesting to note that the acetyl choline receptor 
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displays a dependence on hydrostatic pressure, which becomes 
apparent at 400 bars (44), corresponding to a shift of the tran- 
sition temperature of about 10°K to higher temperatures 1341 . 
This is the approximate difference of the melting temperature 
of biological membranes and physiological temperature. One 
would also expect that an increase in voltage would change the 
state of membranes because each monolayer of the membrane 
possesses a quite high dipole potential of order 200-500mV 
(e.g., B31 ). Antonov et al. Il20l showed that increase in volt- 
age increases the chain melting temperature. Here, we found a 
non-linear current- voltage relationship (Fig. [5]). The increase in 
conductance at higher voltages observed in our BLM measure- 
ments (Fig. [5]l is in agreement with Antonov's findings. How- 
ever, the voltage effect is not well studied. Since the bilayer 
possesses two monolayers with opposing dipole potentials, the 
effect of voltage may be different in the two layers. Impurifica- 
tions and proteins (or their chemical potentials) are further in- 
tensive variables. We have shown previously that fluctuations 
in lipid state are especially enhanced at lipid domain bound- 
aries [46] and at protein interfaces B71 . This makes domain 
boundaries and protein-lipid interfaces particularly susceptible 
for lipid ion channel formation. 

Previously, we have investigated the fluctuation lifetimes of 
lipid vesicles that are identical to the relaxation time scale af- 
ter a small perturbation ll48l I4"9l . We found that relaxation 
time scales are proportional to the heat capacity. Multilamel- 
lar vesicles with very cooperative transitions display relaxation 
times of up to 30 s at the c p maximum. Unilamellar vesicles 
display less cooperative transitions with maximum relaxation 
times around l-2s. If calorimetric profiles are broader (and the 
heat capacity is smaller), one rather expects relaxation times 
in the vicinity of 10-100ms. E.g., it has been estimated that 
the relaxation time scale in the c p maximum of lung surfac- 
tant is about 100ms. This time scale is close to what we find 
here for the channel life times. It has been discussed before 
that channel life times are in fact longer when the heat capac- 
ity is higher |49l . This has been shown in BLMs in a paper 
by Leirer et al. 1381 showing longer mean open times of the 
lipid channels close to the c p -maximum of a DC15PGDMPC 
mixture that correlated with the heat capacity changes. This 
may explain why Antonov et al. [|4) found channel lifetimes on 
the order of seconds at the c p maximum of DPPC membranes. 
As mentioned, this is the relaxation time scale for unilamellar 
DPPC membranes fl49l . Overall, the channel lifetimes seem to 
be related to the fluctuation time scales of lipid membranes. 

Finally, it is tempting to compare lipid channels with protein 
channels. It seems that they display very similar signatures both 
with respect to conductances and to mean channel open times. 
This is striking and raises the question of whether the physics 
behind these events is similar. It is remarkable that Na + chan- 
nels and the acetylcholine receptor are affected by octanol in 
a manner similar to our lipid channels l50l I3TI . The same is 



true for the effect of ethanol that leads to an increase in channel 
life times (data not shown) and which has also been reported to 
have a similar affect on the acetylcholine receptor l52l l53l . In 
the literature, the effect of anesthetics on ion channels has led 
to the notion that anesthesia acts via a specific binding mecha- 
nism to channel proteins, a view that we have opposed in l26l 
on thermodynamic grounds. The similarity between protein 
channel data and lipid membrane conductance deserves to be 
investigated in greater detail in future studies. At present, the 
precise connection between the two phenomena must remain 
an open question. 

Conclusion 

We have shown here that lipid membranes are permeable to 
dyes and ions if one is close to the chain melting regime. The 
microscopic permeation for ions has been shown to be quan- 
tized with conductances and life times similar to those reported 
for proteins channels. Membrane permeation can be influenced 
by changes in intensive thermodynamic variables, in particu- 
lar temperature, voltage and chemical potential of anesthetic 
molecules. Other authors have shown that permeation for ions 
can be influenced by calcium lfl9l [181 . pH changes lfl5l [TBI 
and changes in membrane tension |flT|. The same will be true 
for an increase of hydrostatic pressure. Thus, it must be con- 
cluded that the overall physics of lipid membrane permeation is 
straight forward and easily understood as a consquence of lat- 
eral area fluctuations. Since all biological membranes possess 
lipid membranes (and membrane proteins) it cannot be con- 
clusively excluded that some of the quantized current events 
reported in the literature for proteins are in fact due to lipid 
channels, in particular when considering that biological mem- 
branes exist in a state close to their melting transitions |S]|6]. 
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